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基于 混沌 Duffing 振子 的 BPSK 信号 K-means 聚 类 解 调 方法 


薄 亮 亮 ,， 江 虹 ， 闵 
(西南 科技 大 学 信息 工程 学 院 , 四 川 绵阳 621010) 


摘 要 : 针对 二 进 制 相 移 键 控 (binary phase shift keying，BPSK) 信号 在 低 信 骂 比 下 解 调 误 码 率 较 高 的 问题 ， 提 出 了 一 
种 基于 混沌 Duffing 振子 的 K-means 解 调 方 法 。 该 方法 的 思想 是 根据 混沌 Duffing 振子 系统 对 初 值 的 敏感 性 以 及 对 嗓 
声 的 免疫 特性 ， 低 信 骂 比 下 Duffing 振子 系统 输入 BPSK 信号 时 ,由 于 BPSK 信号 相位 在 0" 和 T 之 间 的 跳 变 从 而 导致 
Duffing 振子 输出 相 轨迹 状态 发 生 改 变 。 针 对 相 轨 迹 状态 的 变化 ,采用 开 -means 聚 类 算法 对 相 轨迹 进行 迭代 求 质 心 ， 根 
据 收敛 后 的 质心 间距 大 小 对 BPSK 信号 进行 判决 解 调 。 仿 真 结果 表明 , 与 现 有 的 几 种 解 调 方法 相 比较 ,基于 混沌 Duffing 
振子 的 BPSK 信号 K-means 聚 类 解 调 方法 在 低 信 骂 比 下 解 调 速 度 、 解 调 精 度 等 方面 都 有 了 较 大 的 提高 。 
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K-means clustering demodulation method for BPSK signal based on chaotic Duffing oscillator 
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Te Abstract: For the binary phase shift keying(BPSK) signal always get high bit error rate(BER) when demodulated in a low 


signal-to-noise ratio(SNR) environment, a K-means clustering demodulation method based on chaos Duffing oscillator has been 


eS proposed to reduce the BER. This study according to two characteristics of Duffing oscillator, not only is it very sensitive to the 


initial value, but also has certain Immunity to noise. The phase of BPSK signal change between 0° and x cause state change in 


站 入 the Duffing oscillator’s phase orbits, and using a K-means clustering method to find the centroid of phase orbit by iteration, then 


> _ BPSK signal can been demodulating by calculate the distance of centroids. Simulation results show that the K-means clustering 


demodulation method has a great improvement in the calculation rate and the demodulation accuracy in a low SNR environment. 
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0 引言 

在 数字 通信 技术 领域 中 ， 由 于 BPSK 调制 信号 具有 良好 的 频带 利用 率 ， 较 强 的 抗 干 扰 性 能 ， 硬 件 易 实现 等 优良 品 性 。 因 此 ， 
BPSK 信和 号 被 广泛 应 用 于 卫星 通信 、 数 字 电视 等 通信 系统 中 。 在 数字 通信 系统 中 ， 由 于 信号 处 于 的 信道 环境 日 趋 复杂 ， 信 道 背 景 
噪声 、 码 间 串 扰 等 因素 的 影响 ， 导 臻 通信 系统 在 接收 端 产生 较 高 的 解 调 误 码 率 。 为 了 提升 接收 端 对 调制 信号 的 接收 解 调 性 能 ， 保 


障 通信 质量 ， 因 此 ， 在 低 信 噪 比 环境 下 降低 信号 解 调 误 码 率 具 有 重要 的 意义 。 
针对 低 信 噪 比 环境 下 BPSK 信号 的 解 调 问题 ， 国 内 外 不 少 学 者 将 非 线性 处 理 方法 引入 到 信号 处 理 领 域 中 ， 文 献 [1.2] 利 用 混沌 
Duffing 振子 对 BPSK 信和 号 进行 处 理 ， 提 出 均 方差 算法 解 调 判决 ， 但 是 对 BPSK 信号 采样 率 过 大 ， 且 没有 考虑 到 BPSK 基带 信号 蔗 
码 间 串扰 ; 文献 [3,4] 基 于 随机 共振 方法 对 2FSK 和 4FSK 信号 进行 解 调 ， 在 没有 考虑 到 基带 信和 号码 间 串扰 的 同时 ， 信 和 号 频率 的 大 
小 也 限制 了 解 调 过 程 。 文 献 [5, 6] 研 究 了 一 种 基于 参数 调整 的 随机 
共振 的 BPSK 信号 解 调 方法 ， 该 方法 在 一 定 程度 上 完成 了 低 信 噪 比 下 BPSK 信号 的 解 调 ， 但 在 求解 过 程 中 算法 的 参数 、 信 号 
频率 以 及 幅度 都 对 解 调 过 程 产生 极 大 的 影响 ， 且 调节 随机 共振 参数 与 信号 参数 使 二 者 匹配 十 分 困难 。 参 考 大 量 文献 发 现 ， 现 有 的 
BPSK 信号 解 调 技术 ， 存 在 以 下 几 个 关键 问题 : a) 现 有 的 一 些 较 先进 的 解 调 方法 在 参数 调节 上 存在 很 大 困难 ; b) 在 解 调 BPSK 信号 
时 ， 大 都 文献 均 忽 略 了 基带 信号 的 码 间 串扰 现象 。 
针对 上 述 存 在 的 问题 ,本文 考虑 了 高 斯 信道 传输 条 件 下 
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信号 码 间 串扰 ,并 设计 符合 信道 传输 的 抗 码 间 症 


扰 的 成 型 滤波 器 。 
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结合 混沌 理论 将 Duffing 振子 应 用 在 BPSK 信号 处 理 中 ， 针 对 混沌 阵子 处 理 后 的 信号 ， 采 用 了 K-means 聚 类 解 调 算法 ， 使 低 信 吕 
比 环境 下 BPSK 信号 解 调 的 误 码 率 得 到 降低 。 


1 ”BPSK 信号 滚 降 调制 及 混沌 Duffing 振子 机 理 简 介 


1.1 BPSK 信号 调制 原理 
BPSK 信和 号 通常 用 相位 0” 和 表示 二 进 制 符 号 信息 “1” 和“-1” 其 时 域 表 达 式 为 


Sosx (DD = Acos(27 J .1+9) (1) 


其 中 : 4 是 信号 的 幅度 ， 大 为 载波 频率 ，g 为 相位 ， 当 发 送 码 元 “-1”, 9 为， 发 送 “1” 时 9 为 0” 。 如 图 1a 所 示 ， 是 未 经 成 型 的 
BPSK 信号 时 域 波形 : 

在 大 多 文献 中 ， 均 是 以 图 1 (a) 中 的 BPSK 信号 来 做 相关 分 析 。 但 是 ， 在 实际 通信 中 ， 由 于 相 邻 码 元 间 会 因为 拖 尾 现象 的 存 
在 而 导致 相互 干扰 (如 图 lb 所 示 )， 故 图 1 (a) 中 的 BPSK 信号 不 能 作为 信息 的 载体 在 信道 中 传输 的 。 因 此 ， 需 要 对 图 1 (a) 中 
的 BPSK 信和 号 进行 成 型 滤波 处 理 ， 达 到 消除 码 间 串 扰 的 目的 。 
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(a) 未 成 型 的 BPSK 信号 时 域 波形 
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(b) 码 间 串 扰 现 象 


(c) 不 同 滚 降 系数 q 对 应 的 冲 
图 1 BPSK 信号 . 码 间 串 扰 、 滚 降 滤 波 器 冲击 响应 时 域 波形 图 
消除 码 间 串扰 常用 的 手段 是 设计 成 型 滤波 器 ， 采 用 具有 低 通 特性 的 冲击 响应 作为 成 型 滤波 器 来 防止 码 间 串扰 串 ， 其 表达 式 如 


名 
二 
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下 : 


2 1 
Sin( 二 cos(QA — COS(QA— 
(7) cos( 元 ) ( 元 ) 


hO= 一 一 一 =Sa() 5 O) 


其 中 : 7 是 采样 时 间 间 隔 ，a 是 升 余弦 深 降 系数 ， 且 满足 0<a<1， 取 wu 值 分 别 为 0，0.5，0.8，1， 得 到 如 图 lc 中 不 同 滚 降 系数 a 
对 应 的 成 型 滤波 器 的 冲击 响应 。 可 见 滚 降 系数 a 越 大 ， 冲 击 响应 函数 h(0 的 拖 尾 衰减 越 快 ， 对 位 定时 的 精度 要 求 越 低 。 成 型 滤波 
器 的 滚 降 系 数 越 大 传输 带宽 增 大 , 带宽 增 大 为 B=(1+a)f, 但 是 带宽 增 大 会 导致 带宽 利用 率 减 小 四 , 带宽 利用 率 为 =R8/B=2/(1+a)。 
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(1) 和 (2) 式 ， 成 型 后 的 BPSK 调制 时 域 表 达 式 为 : 


Somx(D=A(G)ecos(27 大 +OD) (3) 


其 中 ，A(t)= a(t)xh(t)，a(t) 为 码 元 符号 表达 式 。] 
载 频 f=8K Hz， 采 样 频率 f=40f.=320K Hz 深 降 系数 o=0.7， 得 到 如 图 2a 中 BPSK 信和 号 成 型 调制 过 程 ， 图 2b 


BPSK 信号。 


1.2 混沌 Duffing 振子 简介 


混沌 Duffing 振子 具有 初 值 的 敏感 性 以 及 对 噪声 的 免疫 性 。 
另外 ， 混 沌 Duffing 振子 具有 对 噪声 的 免疫 性 ， 即 在 噪声 的 扰动 下 ， 混 六 
低 信 噪 比 环境 下 ， 弱 信号 检测 处 理 和 接收 解 调 带 来 巨大 的 优势 。 本 文 利 ) 


性 完成 低 信 品 比 下 BPSK 信号 的 解 调 。 


Duffing 振子 系统 的 标准 动力 学 行为 方程 
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(a)BPSK 信和 号 成 型 调制 过 程 时 域 波 形 
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(b) 不 同 SNR 下 的 BPSK 信和 号 
图 2 ”BPSK 深 降 调制 过 程 时 域 波形 


对 


述 为 


x+kxtf (xX)=7ycosQxfp) (4) 


其 中 : x 为 Duffing 振子 系统 的 输出 ，k 为 系统 的 阻尼 系数 ，ycos(2mjet) 为 Duffing 振子 的 周期 策动 力 ， 
a 和 4。 均 是 非 线 性 回复 力 参数 M09， 且 满足 a>0,5>0， 故 Duffing 振 


fx) 是 一 非 线 性 回复 力 函 数 ， 
子 方程 表达 式 可 以 为 


其 表达 式 为 fx)=-ax+bx?。 其 中 ， 


x+kx-axtbx? =ycos (2 fn G) 


式 〈6) 亦 可 以 有 如 下 状态 方程 


X=y 


(6) 


y=-B+ar-pe -ycosCzFD 


在 Holmes 型 Duffing 振子 系统 中 ， 取 a=1,b=1， 阻 尼 系 数 广 0.5。 在 不 考虑 载 频 估 计 4 
面值 y 分 别 取 0.8，0.826 和 0.83， 并 用 四 阶 龙 格 - 库 塔 法 解 (6) 式 ， 得 到 Duffing 振子 输出 相 轨 迹 分 别 如 


取 长 度 Len=20 个 的 随机 码 元 ， 基 带 码 元 信号 


混沌 振子 系统 初 值 的 微小 变化 会 引起 系统 输出 
振子 系统 依然 可 以 保持 自 
混沌 Duffing 振子 对 初 值 的 敏感 性 


扁 差 的 情况 下 ， 根 提 


于 混沌 Duffing 振子 的 BPSK 信号 K-means 聚 类 解 调 方法 


名 度 为 0.1V， 码 元 速率 Re=1K baud， 


为 不 同 SNR 下 的 


很 大 的 改变 中; 


~ 


身 丰 富 的 动力 学 特征 0。 这 为 
以 及 对 噪声 的 免疫 特 


7 为 周期 策动 力 的 幅 值 ， 


届 Melnikorv 方法 03]， 


多 


3 吕 


bP (Ca)~ (ec) 所 
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2 duffing r=0.8 混沌 状态 2 duffing r=0.83 周期 状态 
1 1 
0 0 

> > 】 

-1 = -1 NN , 

2 "EE Ox 2 部 0 x 2 

(a)y=0.8 时 的 相 轨迹 (b)y=0.826 时 的 相 轨 迹 (c)y=0.83 时 相 轨 迹 图 


图 3 不 同 策动 力 幅 值 y 对 应 相 轨 迹 图 
3 可 知 ， 随 着 周期 策动 力 信号 幅度 的 增加 ，Duffing 振子 的 相 轨 迹 图 由 混沌 状态 到 临界 混沌 状态 ， 最 后 到 周期 状态 。 因 
此 ， 设 置 合 适 的 周期 策动 力 信号 的 幅 值 《 本 文中 设置 周期 策动 力 幅 值 为 0.826)， 当 加 入 待 测 信号 时 ， 就 可 以 改变 Duffing 振子 输 
出 相 轨 迹 的 状态 。 

在 本 文中 ， 将 待 解 调 的 BPSK 信号 加 入 到 Duffing 振子 系统 中 ， 式 (5) 可 以 写成 


下 


Ey bx =ycos(2xf.1) 0) 
+sn(1) 


Te 


噪声 。 妹 此 ， 


其 中 : sn(t) 是 经 过 传输 信道 引入 噪声 后 的 BPSK 信号 ， 即 sn(t)=Sppsx(t)+n(t)。n(t) 是 均值 为 0, 方差 为 2 的 加 性 高 斯 白 
输入 信 噪 比 可 定义 为 


A-(D) 
| (8) 


bh sk 
SNR, = llogt— loa 


n 


其 中 : Pspst 和 PP; 分别 为 BPSK 信号 功率 和 高 斯 噪声 功率 ， 当 加 入 BPSK 信号 时 ，Duffing 系统 的 方程 可 以 写成 


x+ k x axt+bx =ycos(2xf1) (9) 
+A(f)ecos(2zf .t+09)+n(t) 


ON 对 式 (9) 式 右边 进行 辅助 角 计算 可 得 
> ycos(2z ft) +Alt)cos(2zf.t +9) 


=pP(1)cos(2zxf.t+0(1)) +n(t) (10) 
= BW = Vr +2yA(t)cos p+ A 

en ; 11 
~ 0(1) = arctan Ysin(g) 0 

y+Al(t)cos(9) 
对 于 BPSK 信号 ，92 只 有 0 和 两 种 取 值 ， 故 0(t)=0， 且 : 

_Jz-40D ep=z 
p0-| 和 人 i (12) 


混沌 两 种 状态 的 交 蔡 ， 


cEt 


因此 ， 根 据 周期 策动 力 信号 与 BPSK 信和 号 登 加 后 幅度 p(t) 的 变化 ，Duffing 振子 的 相 轨迹 的 处 于 周期 、 
这 两 种 状态 的 交替 正 由 BPSK 信号 码 元 “1 和 “1” 交 蔡 引 起 的 ， 可 由 图 4 模型 描述 。 


T 
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| 


图 4 Duffing 振子 系统 相 轨 迹 状态 变化 村 


THHI 


图 


根据 图 3 和 4， 系统 输出 相 轨 迹 后 ， 即 可 设计 某 种 判决 机 制 进行 解 调 ， 根 据 混沌 Duffing 振子 输出 


种 K-means 聚 类 解 调 方法 。 


2 BPSK 信号 的 K-means 聚 类 解 调 算法 


K-means 聚 类 算法 属于 简单 的 无 监督 学 习 算法 。 在 K-means 聚 类 中 ， 首 先 ， 将 样 


ChinaX 


聚 类 成 K 个 艇 和 随机 选 ] 


点 ; 然后 ， 计 算 K 个 入 的 质心 点 ， 并 将 计算 得 到 的 K 个 质心 点 进行 重新 聚 类 ;最 后 重复 计算 新 的 聚 类 质心 点 


161。K-means 迭代 求 聚 类 质心 算法 流程 图 如 图 5 所 示 。 


计算 相 轨 迹 上 每 一 个 点 到 
心 的 距离 ， 并 将 该 点 归 
类 到 距 高 近 的 那个 类 中 


vy 


计算 新 的 聚 类 质心 点 


聚 类 质心 点 是 否 改变 从 


图 5 K-means 聚 类 算法 流程 图 

当 Duffing 振子 系统 输入 的 BPSK 信号 不 加 噪声 时 ，1 码 元 〈 相 位 0) 产生 的 相 轨 迹 和 - 
6 a 所 示 ， 当 信 噪 比 SNR=-20dB 时 ，1 码 元 和 -1 码 元 的 相 轨 迹 如 下 图 6b 所 示 。 可 以 看 出 
糙 了 些 ， 并 没有 使 相 轨 迹 状 态 发 生变 化 。 


不 加 噪声 


(不 加 噪声 和 加 噪声 情况 下 , 码 元 1 产生 的 相 轨 迹 图 


人 0 2 和 一 2 


码 元 (相位 x ) 产生 的 相 轨 迹 如 下 图 
噪声 只 是 在 某 种 程度 上 使 相 轨 迹 变 得 粗 
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相 轨 迹 的 特征 ， 本 文 提出 一 


取 K 个 聚 类 质心 
， 直 到 质心 收敛 0 
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2 -2 0 2 
(b) 不 加 噪声 和 加 噪声 情况 下 , 码 元 1 产生 的 相 轨 迹 图 
2， 加 噪声 后 1 码 元 质心 收敛 图 》 加 噪声 后 -1 码 元 质心 收敛 图 
一 相 轨迹 
一 初始 质心 
~ 刀 收 和 敛 后 质心 


多 2 
(C)SNR=-20dB 时 , 码 元 1、-1 码 元 对 应 相 轨 迹 质心 变化 
图 6 码 元 1 和 -1 的 轨迹 及 相 轨 迹 KK-means 质心 收敛 示意 图 
根据 Melnikov 方法 求解 , Duffing 振子 相 轨 迹 图 具有 两 个 鞍点 09, 为 〈 士 10)， 即 不 加 噪声 时 ， 相 轨迹 的 初始 质心 为 ( 士 1.0)， 
对 于 输入 加 噪声 的 信号 , Duffing 振子 系统 输出 的 相 轨 迹 经 过 -means 聚 类 算法 迭代 求 质心 后 , 1 码 元 和 -1 码 元 质心 收敛 后 如 图 6c 
所 示 ， 不 难看 出 1 码 元 对 应 的 相 轨 迹 收 敛 后 的 质心 与 初始 质心 相 比 偏 移 量 较 小 ， 而 -1 码 元 对 应 相 轨 迹 收敛 后 的 质心 偏 移 量 较 大 。 
以 收敛 后 的 两 个 质心 的 水 平 距离 d=K-xi 为 度量 ， 其 中 x1/，xwi 分 别 为 收敛 后 左右 质心 的 横 坐 标 。 在 输入 信 噪 比 为 -134B、-20dB、- 
25dB、-30dB 条 件 下 ， 分 别 取 40 个 样本 1 码 元 和 -1 码 元 ,经 过 K-means 聚 类 算法 后 两 种 码 元 对 应 相 轨 迹 收敛 后 的 两 个 质心 距离 4 
分 布 如 下 图 7 所 示 。 


ye > Gi 0 40 个 1 码 元 | 1.5| : : 全 -40 个 1 码 元 四 
和 ”全 40 个 -1 码 元 | 14 全 40 个 -1 码 元 
“0 10 20 30 40 0 10 20 30 40 
(a)SNR=-15dB 时 ,1 码 元 、-1 码 元 质心 收敛 后 距离 分 布 (b)SNR=-20dB 时 ,1 码 元 、-1 码 元 质心 收敛 后 距离 分 布 
2 - T 2 a 


1.6| 1.6|ee 

1.5| 40 个 1 码 元 1.5 一 一 二 一 40 个 1 码 元 
1.4 | | -全 40 个 -1 码 元 | 1 4| : | -6-40 个 -1 码 元 | 
“0 10 20 30 40 “0 10 20 30 40 


(CJSNR=-25dB 时 ,1 码 元 .-1 码 元 质心 收敛 后 距离 分 布 (JSNR=-30dB 时 ,1 码 元 -1 码 元 质心 收敛 后 距离 分 布 
图 7 不 同 SNR 条 件 下 ，K-means 算法 收敛 后 质心 距离 分 布 图 

从 图 7 不 难看 出 ， 经 过 K-means 聚 类 计算 收敛 后 ，1 码 元 对 应 相 轨 迹 质心 模 坐 标 距 离 范围 约 在 1.76~2.0， 而 -1 码 元 对 应 相 轨 
迹 质心 模 坐 标 距离 范围 约 在 1.48~1.74。 设 置 合适 的 距离 阔 值 ， 利 用 相 轨 迹 数 据 进 行 K-means 聚 类 计算 ， 求 出 收敛 后 的 质心 横 坐 标 
距离 ， 并 与 阔 值 比较 ， 便 能 判别 出 码 元 信息 。 其 步骤 简要 描述 如 下 : 
a) 设 Duffing 振子 输出 相 轨 迹 上 有 nn 个 点 ,设置 初始 质心 为 (+1,0), 计算 相 轨 迹 上 所 有 点 到 质心 (-1,0) 的 距离 di1,q21...di,.….,dn1， 
其 中 da=[Caii+lD2+Oa-0)312， 到 质心 (-1,0) 的 距离 412,422,....,4dz,.….,dn2， 其 中 dio=[(xi2-1)2+Qx2-0)2]12; 

b) 根 据 di 和 diz 的 计算 结果 , 将 距离 (-1,0) 近 的 点 归 为 一 类 , 距离 (1,0) 近 的 点 归 为 一 类 , 这 样 可 将 相 轨 迹 上 的 点 重新 分 成 两 类 ， 
两 类 点 数 分 别 为 nl 和 n2 (nl+n2=n); 

oc) 重新 计算 两 类 点 的 质心 (Xi1,7 了 7)，02, 姓 ) 为 新 的 质心 点 ， 其 中 寅 =(1/n1)Sum(xix2,.xnl)， 了 =(1/n1)Sum(y1y2,...yn1)。 同 理 ， 角 、 
二 亦 可 计算 出 ; 


Wh 
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d) 重 复 9) 中 的 内 容 ， 直 到 质心 点 不 再 偏 移 〈 或 满足 某 个 很 小 的 偏 移 量 ); 


e 计 算 收 敛 后 左右 两 个 质心 的 横 坐 标 之 间 的 距离 =F - zl， 并 


又 ， 使 解 调 端的 计算 量 得 到 降低 。 与 文献 [3] 中 的 随机 共 


与 距离 阔 值 Vw 进行 比较 ， 当 d>Viaq 时 判决 为 1 码 元 ， 反 之 为 - 


数 在 很 大 程度 上 会 影响 最 后 的 计算 结果 ， 并 
的 计算 复杂 度 ， 也 不 需要 算法 参数 与 信号 匹 本 


表 1 


法 、 文 献 [5] 中 提出 的 随机 共振 解 调 法 以 及 传统 的 相 3 
方面 进行 比较 ， 得 到 仿真 结果 如 表格 1 中 所 示 。 


忆 。 


以 上 步骤 分 析 可 知 ，K-means 聚 类 判决 算法 复杂 度 较 低 ， 与 相干 解 调 相 比 ， 在 解 调 端 省 略 了 乘 以 恢复 载波 和 低 通 滤波 等 步 
扑 算 法 相 比 ， 文 献 [5] 中 随机 共振 算法 参数 较 多 ， 且 参数 调节 过 程 复杂 ， 参 
要 求 信 号 自身 参数 与 随机 共振 的 参数 匹配 ， 而 本 文 所 提 方 法 避 开 了 算法 中 参数 带 来 


基于 Matlab 仿真 环境 ， 信 噪 比分 别 设置 为 -234B、-20dB、-15dB、-10dB， 将 本 文 提出 的 基于 Duffing 振子 的 K-means 解 调 方 
F 解 调 法 在 同一 台 计 算 机 上 运行 ,从 仿真 时 间 (算法 复杂 度 )、 解 调 误 码 率 两 个 


K-means 聚 类 解 调 法 在 不 同 SNR 条 件 下 仿真 结 


仿真 参数 设置 : 码 元 长 度 Len=500， 码 元 速率 Ro=1Kbaud 载 频 f=8K Hz， 采 样 频 率 f=40f:=320K Hz， 深 降 系 数 a=0.7， 质 心 距离 闽 值 为 
Vina=1.75。 


信 噪 比 SNR(dB) 


解 调 误 码 率 (%) 


本 文 解 调 法 随机 共振 解 调 法 相干 解 调 法 
-25 16.2 24 30.2 
-20 10 16.2 21.8 
-15 4.4 8 12.6 
-10 1.0 3.8 6 
平均 仿真 时 间 (s) 36.84 58.64 42.76 


BER 为 纵 轴 ， 分 别 计算 本 文 所 提 方 法 、 文 献 [5] 


表 1 不 难得 出 , 在 相同 的 信 噪 比 下 , 本 文 所 提 的 方法 解 调 误 码 率 明 显 低 于 另外 两 种 方法 , 且 计算 速度 也 快 于 另外 两 种 方法 。 
曾 加 仿真 信 噪 比 区 间 ， 每 种 信 噪 比 下 仿真 20 次 ， 其 他 仿真 参数 不 变 。 以 接收 端的 输入 SNR 为 横 轴 ，BPSK 信号 解 调 误 码 率 
的 随机 共振 法 以 及 相干 解 调 法 的 误 码 率 ， 得 到 如 图 8 所 示 的 SNR-BER 曲线 。 由 
图 8 可 知 ， 当 输入 信 噪 比 低 至 -20dB 时 ， 本 文 所 提 解 调 法 误 码 率 约 为 10%， 随 机 ee 16% 和 22%; 


当 输 入 信 噪 比 为 -10dB 时 ， 本 文 所 提 解 调 法 误 码 率 仅 为 1% 左 右 ， 随 机 共振 法 和 相干 解 调 法 误 码 率 分 别 为 4% 和 6%， 通 过 仿真 再 次 
验证 了 本 文 所 提 方 法 在 低 信 噪 比 环境 下 解 调 的 优越 性 。 


35 


-本文 解 调 法 
| i | -e-_ 文 献 [5] 随 机 共振 法 ] 
六 25 i - Se 
由 20H- | | 
加 15| 三 森 
了 
人 
吐 下 es 
0 20 1 0 0 
SNR(dB) 
图 8 K-means 聚 类 解 调 算法 误 码 信 噪 比 - 率 曲线 


3 ”结束 语 


在 数字 通信 中 ， 低 信 噪 比 环境 下 的 BPSK 信号 解 调 十 分 困难 ， 
用 简单 易 行 的 KK-means 聚 类 作为 解 调 算法 。 该 方法 将 
性 迭代 求解 后 ， 得 到 Duffing 振子 的 输出 相 轨 ， 再 经 过 K-means 村 


氏 信 噪 比 下 的 BPSK 信号 作为 混沌 Duffing 振子 的 输入 信号 ， 经 过 


为 了 提升 通信 质量 ， 将 混沌 振子 理论 引入 到 数字 通信 中 ， 并 采 
沌 非 线 


号 让 
到 


案 类 算法 对 相 轨 迹 办 代 求 质心 ， 并 采用 收敛 后 两 质心 


为 判决 依据 ， 进 行 BPSK 基带 码 元 判决 。 仿 真 表明 本 文 所 提 方 法 在 信 噪 比 低 至 -10dB 时 ， 解 调 误 码 率 仅 为 1%， 且 与 其 他 两 种 解 调 


方法 进行 比较 ， 得 出 本 文 所 提 方 法 的 可 行 性 书 


[优越 性 。 
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